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In this article we have studied electrical current and magnetoresistance of the (Fi,Fa) — QD a ~ 
QDb — S system. We have found that the signal of the magnetoresistance can be changed through 
the external potential applied in the ferromagnets. In addition, it is possible to control the current of 
the first ferromagnet (Fi) through the potential applied in the second one (.Fa)- This transistor-like 
behavior can be useful in future technological applications. 



PACS numbers: 73.23Hk, 73.63Kv, 74.45. +c, 74.78Na 

The recent advances in the experimental development 
of nanostructures are mainly aimed at the study of purely 
quantum phenomena and effects based on electron-spin 
properties (spintronics). In the special case of systems 
composed by a ferromagnet and a superconductor it is 
possible to construct spin valves which control the cur- 
rent flow through these systems. Andreev reflection per- 
mits such control, by varying the polarization of the fer- 
romagnet attached to the system, as shown in several 
papers^—. Andreev reflection^ (AR) is a mechanism 
in which a Cooper pair is formed in the superconduc- 
tor from the combination of an incident electron coming 
from the normal metal with energy w and spin a, with 
another electron with energy —uj and spin a. Both elec- 
trons enter the superconductor as a Cooper pair, leaving 
a reflecting hole in the ferromagnetic electrode. Andreev 
states are located within the superconductor gap, where 
no quasi-particles states are available. 

In this work we propose a prototype of a molecular 
transistor by combining two quantum dots with a super- 
conductor and two ferromagnetic electrodes. A schematic 
diagram of the system is shown in figure [T] There are 
two ferromagnetic electrodes, Fi and F2, attached to the 
first quantum dot and a superconductor electrode is con- 
nected to the second one. The dot coupled to the fer- 
romagnetic electrodes (F) is called a, and b is the one 
coupled to the superconductor (S). The superconductor 
has its chemical potential fixed to zero, and independent 
voltage bias are applied to the ferromagnets which are 
called Vi and V2. 




FIG. 1: (Color Online) Schematic diagram showing the 
(F\,F2)-QD\-QD2-S system. The magnetization of Fi is as- 
sumed to be fixed and the magnetization of Fa can be varied 
for an angle 9 with respect to the Fi magnetization. It is ap- 
plied external potentials to F\ and F2 and the superconductor 
lead is grounded. 



The system displayed in Fig. Q] is described by the 
following Hamiltonian: 

H = T-L\ + Ti.2 + T~Ls + Hdqd + 7~Lt, (1) 

where Hi = J2ka( e k ~ — Hi)a\ a a ka is the Hamil- 
tonian of the ferromagnet 1 in which the spin bands 
are shifted by the exchange energy h\. The magneti- 
zation of this ferromagnet is fixed to a direction which 
we call z. The ferromagnet 2 is modeled by H2 — 
J2kA e k - sgn(cr)/i 2 cos6-ni]b\ a b ka - J2 krT h 2 sin 6b\ a b ks 
where 9 is the angle of the F2 magnetization with re- 
spect to z-axis and hn is the exchange energy. The 
superconductor is described by the BCS Hamiltonian: 

Us = E pCT e P s U s P° + E P [ As pt s -;4 + H ; c -]- The c L uan " 
turn dots are considered noninteracting with one level de- 
generated in spin Jidqd = J2a E a ah a <j + J2a E ha h hG . The 
tunneling between the quantum dots and the leads is de- 
scribed by Ht = J2ka MiL'W + H.c] + J2ka ML^W + 

H ' C -] + Ep CT M P <T C ^ + H.C.] + Y.a^abC^Cba + H.C.]. 

In order to calculate the transport properties we have 
used the non-equilibrium Green's function method—. All 
the physical quantities can be cast in terms of the Green's 
function of the dots. By using the Nambu 4x4 notation 
the retarded Green's functions of the quantum dots are 
given by: 

Gr /~irO 1 f-ir + t prO± /~irO fry\ 
aa — ^aa + ^*aa l ab^bb l ab^* a a \ Z ) 
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where we have defined: = ga a (l — 

E^J- 1 and G[° = g&(l " ^M,)^- 

In these equations G r aa is the Green's function of the 
quantum dot a; G£ fc is the Green's function of the quan- 
tum dot b; g^ a and g£ b are the Green's functions of the 
dots a and b isolated from the electrodes; t a j, describes 
the coupling between the dots; H r F = E£ + H r 2 and H r s 
are the retarded self-energies which describe the coupling 
of the dots with the ferromagnetic and superconductor 
electrodes, respectively. Explicitly these self-energies are 
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written as, 
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with A a = r lCT + c 2 r 2ff + s 2 r 2ff , b = sc (r 2T - r 2; ), 

s = sin 0/2 and c = cos 9/2. We also have defined r j(T = 
27r|^| 2 A^ iCT , (with i = 1,2) as the coupling strength, with 

being the tunneling amplitude and Ni a the density of 
states for the ferromagnet spin a band. 

The retarded/ advanced self-energy of the supercon- 
ductor is given by, 



in which the Fermi functions are defined as fi = f(ui — 
eVi) and fi = (uJ + eVi) for i = 1,2 and fi = f(uj) if 
i = s. The electrical current can be calculated as the 
number variation in any point of the circuit. Here we 
choose the variation at the ferromagnet electrodes such 

as Ii = —e ^^E 1 ) w ith i = 1, 2. The total current is given 
by the sum of the current from the electrode Fx and the 
electrode F 2 , i.e., I = I\ + J 2 . The final result is: 



I = - I MMiih - A) + A 12 (fx - / a ) 
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, with the imaginary part accounting for An- A 22 = (r 2t c 2 + s 2 T 2l ) 2 \G r aaM \ 2 +c 2 s 2 {G r aaM [G r aa ^ 12 Y 



where T s = 27r|t s | 2 A^ s , with N s being the density of states 
of the superconductor in the normal state and p is the 

modified BCS density of states p(uj) - 1 ; 

uj9(A-\lj\) 

WA 2 - lo 2 
dreev states within the gapi£. 

Besides the retarded and advanced Green's functions, 
it is necessary to obtain the "lesser" Green's functions 
for the quantum dots, which are obtained through the 
Keldysh equation. By using the equation of motion tech- 
nique we can obtain an exact self-energy expression which 
allows us to write for the dot a: 

G< (w) = G^»£<»GLM, (6) 
with the "lesser" self-energy S^ a : 

£<» = £<H + ti b Gl° b X<(Lo)Gf b (cu)t ab . (7) 

Correspondingly, the Keldysh equation for quantum 
dot b is given by: 
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with the "lesser" self-energy S^ b : 

S< b ( w ) = £<H + ti 6 G^S<(c)G^( w )t ab . (9) 

In the equations (JSJ) through © the "lesser" self- 
energies Tip — Tf + T 2 and are obtained by the 
fluctuation-dissipation theorem Hf = Fi(uj)[Ef — S[], 
where i = 1, 2 or s. The Fermi matrix is given by, 



A 21 = r u [(r 2t c 2 + s 2 r 2i )|GL, 12 | 2 
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The current is composed by four contributions which 
can be interpreted as discussed by the reference HO- The 
terms A\\{f\ — fi) and A 22 (f 2 — f 2 ) represent the direct 
AR in the electrodes F\ and F 2 , respectively; A\ 2 {fx — fi) 
represents the crossed Andreev reflection (CAR) of an 
electron of F\ as a hole of F 2 . The term A 2 \(f 2 — fx) 
represents a similar process for an electron of F 2 as a 
hole of Fx. 

The magnetoresistance of the system is defined as, 
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FIG. 2: (Color Online) ARMR, I P and Iap curves, (a), (b) and (c): V x = F 2 = 0.50; (d), (e) and (f): Ti = 0.20 and T 2 = 0.80. 
Fixed parameters: T s = 0.30, t a t ~ 0.20, Pi = P2 = 0.95, fcjT = 0.01. All the parameters are expressed in superconductor 
gap units. 



in which I A p = 1(9 = tt) and I P = 1(9 = 0). 

The definition (fT2j) is different from the usual since we 
use the absolute value of the currents. This definition 
allows us to compare the amplitude of the currents in 
terms of the bias of each electrode. In this system the 
current is conducted by crossed Andreev reflections, as 
a result, the sign of current in each ferromagnetic ter- 
minal is linked to the averaged chemical potential of the 
two leads. Thus, it contains the case that V\ > and 
V2 < but I > 0. Therefore, by using the definition 
(|12p we can define which current is larger through the 
signal of ARMR even in the cases where we consider the 
dependence of ARMR with the bias V\ or V2. 

In the figure [21 we consider the dependence of the 
ARMR with V\ for different values of V2. The corre- 
sponding Ip and Iap curves are also shown. We analyze 
the influence of the symmetry of the coupling with the 
ferromagnets. 

It can be observed from the figures ^Bp and &p that Ti 
limits the amplitude of the current Ip and T2 shifts Ip 
on the current axis. On the other hand, by comparing 
the figures [2): and [2f, we note that the current Iap is 
almost independent on the relation between Ti and IV 
This is explained by the value of the polarization, which 
is equal to 0.95 for both electrodes. In this case the total 
current is almost carried by the CAR, which picks up one 
spin-up electron from F\ and another spin-down electron 
from F2- Since the total current must be unpolarized by 
the superconductor, it is limited by the electrode with 



low injection of electrons. Therefore, the global behav- 
ior of the magnetoresistance follows the variations of the 
current Ip. 

From the definition (TT2"j) , the current Iap > Ip corre- 
sponds to positive values of ARMR. When the polariza- 
tion is close to unity, the usual situation is to find positive 
values of ARMR since the current Iap is carried by the 
CAR besides the direct AR, which is the only channel 
available to carry the current Ip. In fact, as shown in 
figure Hk, the ARMR curve for Ti = T2 is positive for 
all values of the bias. For the curves with V2 ^ there 
are two peaks, localized at V\ f=s ±0.22, corresponding to 
the value where Iap is close to zero. 

In figure is shown the ARMR curves for T% > T\. 
When V2 = the ARMR is almost constant for all values 
of Vi . For V 2 > the ARMR displays a step-like behav- 
ior with positive values for V± > and negative values for 
V\ < 0. If we change the sign of V% we have the opposite 
case, i.e, ARMR > for Vi < and ARMR < for 
Vi > 0. The results shown in figure [5] indicate that one 
can control the sign of the magnetoresistance of the sys- 
tem through external parameters (Vi and V2) when the 
coupling of F\ and F% with the dot a arc different from 
each other. 

The curves for Iap, shown in figures and[2f, display 
some interesting features. Unlike the current Ip, the shift 
of Iap along the current axis is related with the bias 
rather than the coupling constants Ti and T2- In fact, 
the high polarization values P\ = P2 = 0.95 reduce the 
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FIG. 3: (Color Online) Current through the terminal Fi for 
V2 = —0.30 (red curve) and V2 ~ 0.30 (black curve). Fixed 
parameters: Fi = 0.20, F 2 = 0.80, T s = 0.80, t ab = 0.20, 
Pi = P2 = 0.95, V 2 = 0.30, k B T = 0.01. All the parameters 
are expressed in superconductor gap units. 

contribution from the direct AR and the current Iap is 
dominated by the CAR process. In order to clarify this 
point we consider the current formula (ITT]) in the case 
of zero temperature, where the Fermi functions become 
step functions. If we also neglect the direct AR, I{V\, V 2 ) 
can be written as, 

e f +Vl 

I{Vi,V 2 ) = - [A 12 {u) + A 21 {-u)] dta (13) 
h J-v 2 
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